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Abstract. In an accompanying paper [J. A. Sauls and M. Eschrig, Vortices in 
chiral, spin-triplet superconductors and superfluids, arXiv:0903.0011J we have studied 
the equilibrium properties of vortices in a chiral quasi-twodimensional triplet super- 
fluid/superconductor. Here we extend our studies to include the dynamical response 
of a vortex core in a chiral triplet superconductor to an external a.c. electromagnetic 
field. We consider in particular the response of a doubly quantized vortex with a 
homogeneous core in the time-reversed phase. The external frequencies are assumed 
to be comparable in magnitude to the superconducting gap frequency, such that the 
vortex motion is non-stationary but can be treated by linear response theory. We 
include broadening of the vortex core bound states due to impurity scattering and 
consider the intermediate clean regime, with a broadening comparable to or larger than 
the quantized energy level spacing. The response of the order parameter, impurity self 
energy, induced fields and currents are obtained by a self-consistent calculation of the 
distribution functions and the excitation spectrum. Using these results we obtain the 
self-consistent dynamically induced charge distribution in the vicinity of the core. This 
charge density is related to the nonequilibrium response of the bound states and order 
parameter collective mode, and dominates the electromagnetic response of the vortex 
core. 
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1. Introduction 



Transport and optical properties in type II superconductors are closely related to the 
response of Abrikosov vortices [Ij to external electromagnetic fields. Abrikosov vortices 
that form a lattice in sufficiently strong magnetic fields give rise to a finite resistance 
when a voltage is applied to the superconducting sample. This fiux fiow resistance 
appears, once collective vortex pinning [2] is overcome, as a result of motion of the 
vortex lattice in an electric field [3] with a resistance independent on the field strength 
for small fields [H O [6] . In superconductors with moderate impurity disorder the vortex 
lattice moves in a constant electric field with constant velocity i^l predominantly in the 
direction transverse to the transport current, 

Vl = ^Jtr X ^0 (1) 

with #0 = {hc/2\e\)nB being the magnetic fiux quantum multiplied with the unit 
vector Ub = B.^^/B^^ that describes the direction of the averaged over the vortex unit 
cell magnetic field , -Bav, and /i quantifies the vortex mobility. This viscous fiux fiow 
is in contrast to superclean superconductors or to neutral superfiuids, where vortices 
move predominantly together with the superfiuid velocity. The electric field averaged 
over the unit cell, ^^av is given by 

E^^ = -Sav X -Ul. (2) 

c 

On a phenomenological level, vortex motion has been discussed within the Bardeen- 
Stephen model [7|, the Nozieres-Vinen model [S], and within time-dependent Ginzburg- 
Landau theory [9]. The microscopic theory for fiux fiow resistance has been pioneered by 
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Gor'kov and Kopnin [HE] and by Larkin and Ovchinnikov [H El [101 HH 112] . In extremely 
clean superconductors vortices move predominantly in direction of the transport current, 
giving rise to an electric field perpendicular to the transport current, and thus a Hall 
effect [l3l[Tll[l5l[T6l[I71[T8l[T9l[2n]. 

Vortices also provide valuable information about the nature of low lying excitations 
in the superconducting state. In clean s-wave BCS superconductors the low-lying 
excitations in the core are the bound states of Caroli, de Gennes and Matricon |21j . 
These excitations have superconducting as well as normal metallic properties. For 
example, these states are the source of circulating supercurrents in the equilibrium 
vortex core, and they are strongly coupled to the condensate by Andreev scattering 
[22| |23]. Furthermore, the response of the vortex core states to an electromagnetic field 
is generally very different from that of normal electrons. 

In the following we will discuss the response of a vortex to an a. c. electromagnetic 
field with frequency uj. We are interested in frequencies u that are comparable to the 
frequencies set by the superconducting gap, A/h, i.e. 1 GHz — 1 THz. In this region 
the vortex motion is non-adiabatic and is not described by the low-frequency limit of 
flux flow motion. We consider layered superconductors with a weak Josephson coupling 
between the conducting planes. In this case the vortices are quasi-two-dimensional 
objects, called pancake vortices. The peculiarities of the vortex phases in layered 
superconductors has been reviewed in Ref. [21] (for recent developments see also [25]). 

Disorder plays a central role in the dissipative dynamics of the mixed state of 
type II superconductors. Impurities and defects are a source of scattering that limits 
the mean free path, i, of carriers, thus increasing the resistivity. Defects also provide 
'pinning sites' that inhibit vortex motion and suppresses the flux-flow resistivity [2]. 
However, for a. c. fields even pinned vortices are sources for dissipation. The magnitude 
and frequency dependence of this dissipation depends on the electronic structure and 
dynamics of the core states of the pinned vortex. Experimental studies of dissipation in 
vortex cores include Refs. [26l (23 128] . 

When discussing vortex cores states in the context of impurity disorder, one needs 
to distinguish three degrees of disorder, which we summarize in Table [1] In the dirty 
limit, h/r ^ A, the the Bardeen-Stephen model [7] of a normal-metal spectrum 
with the local Drude conductivity in the core provides a reasonable description of the 
dissipative dynamics of the vortex core. The opposite extreme is the "superclean limit" , 

Table 1. Three degrees of disorder arc distinguished for vortices in superconductors 
with impurities. Here, l/r = Vf/i, with the normal state Fermi velocity Vf. 



superclean case 


moderately clean 


dirty limit 










— « - « 

Ef T 


T 



Charge Dynamics of Vortex Cores in Layered Chiral Triplet Superconductors 4 

h/r <^ / Ef (where Ef is the Fermi energy), in which the quantization of the vortex- 
core bound states must be taken into account. In this hmit a single impurity and its 
interaction with the vortex core states must be considered. The a.c. electromagnetic 
response is then controlled by selection rules governing transition matrix elements for 
the quantized core levels and the level structure of the core states in the presence of an 
impurity |29l [30l [311 |32l [331 [31]. In the case of d-wave superconductors in the superclean 
limit, minigap nodes in the spectrum of bound states lead to a finite dissipation from 
Landau damping for T ^ [35] . 

The superclean limit is difficult to achieve even for short coherence length 
superconductors; weak disorder broadens the vortex core levels into a quasi-continuum. 
We investigate the intermediate-clean regime, A^/Ef <^ h/r <^ A, where the discrete 
level structure of the vortex-core states is broadened and the selection rules are broken 
due to strong overlap between the bound state wave functions. However, the vortex 
core states remain well defined on the scale of the superconducting gap, A. In this 
regime we can take advantage of the power of the quasiclassical theory of nonequilibrium 
superconductivity [361 1311 [Ml [391 [10] ■ Previous results for s-wave and d-wave vortices in 
layered superconductors [lOl [HI [121 133] have shown that in the the intermediate-clean 
regime electrodynamics of the vortex state is nonlocal and largely determined by the 
response of the vortex-core states. Transitions involving the vortex-core states, and their 
coupling to the collective motion of the condensate requires dynamically self-consistent 
calculations of the order parameter, self energies, induced fields, excitation spectra 
and distribution functions. In particular, it has been found that branches of localized 
vortex bound states that cross the chemical potential (so-called anomalous branches) 
are of crucial importance for the dissipation in the vortex core [H]. The relaxation of 
localized excitations takes place via their interaction with impurities. The importance 
of anomalous branches for the fiux fiow resistance in moderately clean superconductors 
has been clarified by Kopnin and Lopatin |17j . 

An interesting aspect of vortex dynamics is the charge dynamics associated with 
moving vortices. It is known that even for a static vortex charge transfer takes place from 
the vortex core region to outer regions of the vortex, leading to a depletion of charge 
carriers in the vortex core regions [I8l[43[45l[46l[47l[48l[4g[50l[5ll[5a[5a 
Associated with this charge transfer is a spatially varying electrostatic potential that 
is typically of the order of A jEj. The smallness of the parameter AjEj means that 
the electrostatic potential can be neglected in the self-consistent determination of the 
spatially varying order parameter and magnetic field. 

There are two main sources for vortex charges. The first contribution comes from 
the lowering of free energy by the pairing, which leads to a force on the unpaired electrons 
in the vortex core by the condensate around the vortices [56] • This force varies on the 
coherence length scale. It arises due to the change of the chemical potential as function 
of the modulus of the order parameter, and is in general proportional to d'yjdn and 
dTc/dn, where 7 is the normal state specific heat coefficient, Tc the superconducting 
transition temperature, and n the density of conduction charge carriers. It includes 
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contributions from the condensation energy -Econ, which is proportional to 'jT^- For 
example, in a Gorter-Casimir two-fluid picture [57], using an extended version [58] of 
the Ginzburg-Landau approximation [59] the corresponding contribution reads [51] 



— ^ + — yj;; (3) 

where ns{r)/n is the condensate fraction of the density. This term is referred to 
as the thermodynamic contribution to the electrostatic potential [GQl [61]. A second 
contribution comes from the inertial forces and the Lorentz force that act on the 
circulating supercurrent around the vortex [621 [63|, |6l] . This force decays on the scale of 
the London penetration depth, and is comparable in magnitude with the thermodynamic 
contribution near the vortex core. It arises due to the kinetic energy of the superfiow, 
and also leads to a depletion of charge carriers in the vortex core region. In the inter- 
vortex regions, where the supercurrent momentum is determined by the phase gradient 
of the order parameter, Vt9, the corresponding contribution to the electrostatic potential 
reads [M] 

e^Bernouiii(r) ^ 77^ {hV^{r) - e*A{r)f (4) 
n 2m* 

where e* = 2e and m* are the effective charge and mass of the Cooper pair, and A 
is the vector potential. The corresponding electrostatic potential is referred to as the 
Bernoulli potential. The electrostatic potential resulting from both contributions leads 
to a charge depletion in the vortex core of order 

5<5static ~ e ^-^^ . (5) 

In general, the electrostatic potential is determined by a screened Poisson equation 
[65]. However, the spatial variations of n{r), A{r), and A{r) occur on long-wavelength 
superconducting scales, while screening takes place on the short-wavelength Thomas- 
Fermi screening length, Atf- Thus, to leading order in the quasiclassical expansion 
parameter, Xtf/^q, the superconductor maintains local charge neutrality [HI [66]. The 
first corrections to the charge density are then given by 

- V2^ = 47rp(r). (6) 

with p(r) = n{r) — hq, and hq is the charge density in the normal state. We note that 
a possible charge pileup in a region of size A|.p around the vortex center induced by the 
presence of the superconducting state is expected to be small, because the corresponding 
overlap region between such a charge distribution and the superconducting order 
parameter is small by the ratio (Atf/^o)^ and in addition the order parameter is small by 
a factor Axp/^olil Consequently, any possible charge induced within the area A^p in the 
vortex center must be small by at least the ratio (Atf/Co)^- This allows us to neglect such 

I This last statement needs to be modified if vortex core slirinking takes place at low temperatures in 
superclean materials, however the general conclusions will still hold. 
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charges and leads to an effective decoupling between charge variations on the Thomas- 
Fermi screening length and the superconducting order parameter. Charge ffuctuation 
on the Thomas-Fermi length scale are entirely determined by the normal state charge 
ffuctuations, and the leading order charge variations induced by the superconducting 
state are varying on the long-wavelength (superconducting) scales. 

The entire picture discussed so far for the equilibrium vortex lattice changes 
dramatically when considering a time-dependent perturbation on a vortex. In order 
to reduce the Coulomb energy associated with the charge accumulation an internal 
electro-chemical potential, ^{r;t), develops in response to an external electric field [6]. 
This potential produces an internal electric field, £''"'(r;t), which is the same order of 
magnitude as the external field. Even though the external field may vary on a scale 
that is large compared to the coherence length, ^o, the internal field develops on the 
coherence length scale. The source of the internal field is a charge density p(r, t) that 
accumulates inhomogeneously over length scales of order the coherence length. It is 
necessary to calculate the induced potential self consistently from the spatially varying 
order parameter, spectral function and distribution function for the electronic states in 
the vicinity of the vortex core. An order of magnitude estimate shows that to produce 
an induced field of the order of the external field, the dynamically induced charge is of 
order 

5<5 dynamic ~ 6 — (7) 

where 6v^ ~ cE'^'^^^qA/u! is the typical energy scale set by the strength of the external 
field. This charge density accumulates predominantly in the vortex core region and 
creates a dipolar field around the vortex core [3, HI] . For a pinned vortex the charge 
accumulates near the interface separating the metallic inclusion from the superconductor 
|43] . The charge in Eq. ([7]) exceeds the dynamical charge that would result from a 
dynamic motion of the vortex charges in Eq. ([5]), provided 6v^ ^ A^/Ef. In this 
case, we can neglect the contributions coming from the response of the charges that are 
present already for a static vortex, and concentrate on the new contributions Eq. ([7]). 
Note that the dynamical charges are absent in the superclean case, where the electric 
field is purely inductional [T6] . 

In the next two sections we provide discussion of the expansion in small 
quasiclassical parameters, and a summary of the nonequilibrium quasiclassical 
equations, including the transport equations for the quasiparticle distribution and 
spectral functions, constitutive equations for the order parameter, impurity self-energy 
and electromagnetic potentials. We then apply this theory to the response of a vortex 
in a chiral, spin-triplet superconductor to an a.c. electromagnetic perturbation. The 
equilibrium properties of such a vortex have been discussed in detail in an accompanying 
paper, Ref. [67]. In particular, we discuss the a.c. response of singly quantized and 
doubly quantized vortices in Sec. HI 
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2. Expansion in small parameters 

The physics of inhomogeneous metals and superconductors described by the 
quasiclassical approximation is governed by well defined small expansion parameters. 
We will assume that there is one such parameter that describes all small quantities 
in the system, and will assign to this parameter the notation small [HHl [Ml EOl ITT] . 
The typical microscopic length scales of the problem, in short denoted by oq, are the 
Bohr radius as, the Fermi wave length A/, and the Thomas- Fermi wave length Atf- 
The mesoscopic, superconducting length scales are the coherence length ^q, and the 
penetration depth A. Correspondingly, large energy scales are the Fermi energy Ef, and 
the Coulomb energy e^/ao, whereas small, quasiclassical energy scales are the gap. A, 
the energy of the external perturbation, hu, and the energy related to the transition 
temperature, /crT'c. We assume that spatial variations near the vortex core are on the 
scale .^0; aiid time variations on the scale I/uj. The normal state density of states at 
the Fermi level, Nf, is of order Nf ~ l/Efa^. With this notation we have the following 
assignments that we need to estimate the electromagnetic fields, charges and currents: 



smalP 


small^ 


smalP 


e^Nf ~ \ 
al 






hvf ~ 


Vf ao 
c ^ \ 






hvf ■'V ^ A 


SO 




hdt ~ huj 


1 ao hid 
-Of ~ 



2.1. Electrostatic Fields and Potentials 

In equilibrium, the electrostatic fields are given in terms of the potentials hy E = — V<P, 
S = V X A, and fulfill the Maxwell's equations V-S = 0, VxB = ^j, V ■E = inp, 
^ X E = 0. The sources for the static fields in the superconductor are the Meissner 
currents j and the small charges p discussed in the introduction. The Meissner current 
density in a superconductor can be estimated by noting that ~ e^VfNfA/c ~ 
Z\/aoA ~ smalP. The electrostatic potentials induced in the superconducting state are 
of order e<P ~ ^^^f ~ ^^/-^/^o; ^ind estimating Ef ~ VfPf/2 ~ hvf/Xj ~ A^o/ao we 
arrive at ~ aoA/^Q ~ smalP. The averaged static magnetic field Bq is given by 
the condition that one fiux quantum penetrates the vortex unit cell. In contrast, the 
variation of the magnetic field on the coherence length scale, SBq, is a factor C,o/X smaller 
than the averaged field. This leads to the following estimates for physical parameters: 
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In particular, for the free energy density the contributions from the Coulomb energy 
and the magnetic field energy are 

^J^r^A^NSr^small', A^Nf^ smalf . (8) 

2 A 2 4o 

For the magnetic field energy the relevant quantity considered here is the contribution 

that results from the deviations 6Bq from the averaged field Bq and its corresponding 

current Sj = -^'V x SBq. Note that in the strong type II limit and for short coherence 

length superconductors (this is the case for example for cuprate superconductors, where 

^o/A ~ 0.01 and ao/^o ~ 0.2), the scales can conspire such that (Co/A)^ is of the 

same order as (ao/^o)'*- In this case, the Coulomb energy contribution might become 

important for the vortex lattice structure [72]. 



2.2. Electrodynamical Response 

For the dynamical response it is useful to split the vector potential into transverse and 
longitudinal parts, 

A = Al + At, V • At = 0, V x Al = 0. (9) 

The longitudinal part can be written as Al = ^VC, and can be gauged away by 
changing the phase of the order parameter 6 into 6 + ( and the electrical potential <P 
into (p' = <P + ^dtC- Thus, we use a gauge where A' = At. The fields are given in 
terms of the potentials as 

B = VxA', E = --dtA' - V <P' , (10) 

c 

and the Maxwell equations read (with the corresponding estimates of the various terms) 
-V'^>' = 47rp, -^^A' + \d^A' + -dtV^' = —j . (11) 

smalfi smalP smaW^ smaW^ 

We drive the superconducting vortices out of equilibrium with an external a.c. 
electromagnetic field, SA^'^^it) = ^AoC"*"^*, of frequency u ~ A/h. The perturbing 
potential then is of the form 6v^^ = —^Vf- 6A^^{t). For the linear response 
approximation to hold, we assume that 6v^^ ~ hooS with 5 a small parameter that defines 
the region of applicability of linear response. This means that for small frequencies the 
vortex oscillation amplitude will be ~ ^o^; and the vortex velocity ~ Vfifiu/ A)5. We 
can then estimate the fields to be of order: 
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Below, wc use these estimates to simplify the coupled system of Maxwell's equations 
and transport equations for the superconducting charge carriers. In particular, we note 
that because 

V ■ i ~ ■ o ~ small^ ■ S, dtp ~ . .„ ~ smaW^ ■ 5, (12) 

the condition of local charge neutrality is fulfilled in linear response to within two orders 
in small. Furthermore, the terms 

\dlA + -dtV^' - small^ (13) 
c^ c 

can be neglected compared to the terms V^A' ~ small^ and —j ~ small^. 



3. Nonequilibrium Transport Equations 

The quasiclassical theory describes equilibrium and nonequilibrium properties of 
superconductors on length scales that are large compared to microscopic scales (i.e. 
the lattice constant, Fermi wavelength, kj^, Thomas-Fermi screening length, Atf, etc.) 
and energies that are small compared to the atomic scales (e.g. Fermi energy, Ef, plasma 
frequency, conduction band width, etc.). Thus, the expansion parameter small defines 
the limits of validity of the quasiclassical theory. In particular, we require kf^Q ^ 1, 
ksTc/Ef <^ 1 and hu <^ Ef, where the a.c. frequencies of interest are typically of order 
A/h ^ k^Tc/h, or smaller, and the length scales of interest are of order the coherence 
length, ^0 = hvf /2nkBTc, or longer. Hereafter we use units in which h — kB = I, and 
adopt the sign convention e = — |e| for the electron charge. 

In quasiclassical theory quasiparticle wavepackets move along nearly straight, 
classical trajectories at the Fermi velocity. The classical dynamics of the quasiparticle 
excitations is governed by semi-classical transport equations for their phase-space 
distribution function. The quantum mechanical degrees of freedom are the "spin" and 
"particle-hole degree of freedom", described by 4x4 density matrices (Nambu matrices). 
The quantum dynamics is coupled to the classical dynamics of the quasiparticles in phase 
space through the matrix structure of the quasiclassical transport equations. 
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The nonequilibrium quasiclassical transport equations [SHI EH EH EE [ID] are 
formulated in terms of a quasiclassical Nambu-Keldysh propagator g{pj^, R; e,t), which 
is a matrix in the combined Nambu-Keldysh space [73], and is a function of position 
R, time t, energy e, and momenta on the Fermi surface. We denote Nambu-Keldysh 
matrices by a "check" , and their 4x4 Nambu submatrices by a "hat" . Thus, the Nambu- 
Keldysh matrices for the quasiclassical propagator and self-energy have the form, 

'=[o g^j ' (14) 

where ^^'^'^ are the retarded (R), advanced (A) and Keldysh (K) quasiclassical 
propagators, and similarly for the self-energy functions. Each of these components 
of g and a are 4x4 Nambu matrices in combined particle-hole-spin space. For a review 
of the methods and an introduction to the notation we refer to Refs. [I2l EH |69] . In 
the compact Nambu-Keldysh notation the transport equations and the normalization 
conditions read 

[{e+^Vf A)f3-eZo^l-A^i-i>^i-ai,g]^ + iVfVg =0, (15) 

g(^g = -nn , (16) 

where the commutator is [A, B](^ = A ^ B — B ^ A, 

A ® 5(e, t) = ei(^^"«f -^*"«?)i(e, t)5(e, t) . (17) 

The vector potential, A{R;t), includes Ao(-R) which generates the static magnetic 
field, Bq{R) = V X Aq{R), as well as the non- stationary vector potential describing 
the time-varying electromagnetic field; A^{{pj:, R; t) is the mean-field order parameter 
matrix, z/^f (pj, -R; t) describes diagonal mean fields due to quasiparticle interactions 
(Landau interactions), and (5"i(p/, R, e, t) is the impurity self-energy. The electrochemical 
potential ^{R;t) includes the field generated by the induced charge density, p{R;t). 
The coupling of quasiparticles to the external potential involves virtual high-energy 
processes, which result from polarization of the non-quasiparticle background. The 
interaction of quasiparticles with both the external potential ^ and the polarized 
background can be described by coupling to an effective potential Zq^ [6H]. The 
high-energy renormalization factor Zq is defined below in Eq. (1251) . The coupling 
of the quasiparticle current to the vector potential in Eq. (fTSl) is given in terms of 
the quasiparticle Fermi velocity. No additional renormalization is needed to account 
for the effective coupling of the charge current to the vector potential because the 
renormalization by the non-quasiparticle background is accounted for by the effective 
potentials that determine the band structure, and therefore the quasiparticle Fermi 
velocity. 

In quasiclassical theory the description in terms of the variables (pj, R; e) is related 
to the (p, R) phase-space description by the transformation of distribution functions 
f{Pf, R, e, t) to n(p, R, t) = f{pf, R, e(p, R, t),t), with e(p, R, t) = Vf{pj) ■ (p - p^) + 
u^iiPf, R,t) + eZo<P{R,t) - -^Vf{pf) ■ A{R,t), see Ref. [68]. 
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3.1. Constitutive Equations 



Equations ( 11511161) must be supplemented by Maxwell's equations for the electromagnetic 
potentials, and by self-consistency equations for the order parameter and the impurity 
self-energy. We use the weak-coupling gap equation to describe the superconducting 
state, including unconventional pairing states. The mean field self energies are then 
given by, 

A^^iPf,R;t) = Nf £"|i(V(p^,p'^)/^(p'^,i?;e,t)) , (18) 

u^^{Pf,R;t)=Nf £"^(A(p^,p'^)r(p},i?;6,t)) , (19) 

A l,{pj, R;t) = , z>^f (p^, R,t)=0 . (20) 

The impurity self-energy, 

ai{pj,R,e,t) = rii i{pj,pj-, R;e,t) , (21) 

is specified by the impurity concentration, n^, and impurity scattering t-matrix, which 
is obtained from the the self-consistent solution of the t-matrix equations, 

i{Pf,p'f, R, e,t) = u{Pf,p"f) + Nf{u{pf,p'j) ® g{p'f,R, e,t) ® i{p'f,p},R, e,t)) , (22) 

where the Fermi surface average is defined by 

and Nf is the average density of states on the Fermi surface. The Nambu matrix 
jK (^g^) is ^\iQ off-diagonal (diagonal) part of in particle-hole space. The other 
material parameters that enter the self-consistency equations are the dimensionless 
pairing interaction, NfV{pf,p'j-), the dimensionless Landau interaction, NfA{pj,p'j-), 
the impurity concentration, ni, the impurity potential, u{ppp'jr), and the Fermi surface 
data: (Fermi surface), Vf(j)j) (Fermi velocity). We eliminate both the magnitude 
of the pairing interaction and the cut-off, in favor of the transition temperature, Tc, 
using the linearized, equilibrium form of the mean- field gap equation (Eq. (ITSl) ). 

The quasiclassical equations are supplemented by constitutive equations for the 
charge density, the current density and the induced electromagnetic potentials. The 
formal result for the nonequilibrium charge density, to linear order in A/Ef, is given in 
terms of the Keldysh propagator by 

p^'\R-t)=eNf £'^ii(Z(p'^)Tr [^^(p}, e, t)] ) ~ 2e'NfZo<I>{R;t) , (24) 

with the renormalization factors given by 

Z{pj) = l-{NjAip'f,p^)), Z, = {Z{pf)). (25) 

The high-energy renormalization factor is related to an average of the scattering 
amplitude on the Fermi surface by a Ward identity that follows from the conservation 
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law for charge [HS]- The charge current induced by A{R;t), calculated to leading order 
in A/Ef, is also obtained from the Keldysh propagator, 

= eNf J ^Tr{vf{p'^)ng''{p'j, R;e,t)) . (26) 

There is no additional high-energy renormalization of the coupling to the vector 
potential because the quasiparticle Fermi velocity already includes the high-energy 
renormalization of the charge-current coupling in Eq. [261 Furthermore, the self- 
consistent solution of the quasiclassical equations for ensures the continuity equation 
for charge conservation, 

dtp^'\R;t) + V-j^'\R;t) = 0, (27) 

is satisfied. An estimate of the contribution to the charge density from the integral in 
Eq. (HM leads to the condition of "local charge neutrality" [6l [66]. A charge density 
given by the elementary charge times the number of states within an energy interval 
A around the Fermi surface implies 2eA^/Zi. Such a charge density cannot be 

maintained within a coherence volume because of the cost in Coulomb energy. This 
can be seen in the estimate in Eq. (fTTI) . where up to third order in the parameter 
small the charge fluctuations p are suppressed. This suppression, and the associated 
suppression in Coulomb energy, must be taken care of by requiring the leading order 
contribution to the charge density to vanish: i.e. p^^\R]t) = 0. Thus, the spatially 
varying renormalized electro-chemical potential, Zq^, is determined by 

2eZoHR,t) = |_^'^^Tr(Z(p})^^(p'^,i?;e,t)). (28) 

The continuity equation implies V ■ j^^\R; t) = 0. We discuss violations of the charge 
neutrality condition f[28|) . which are higher order in A/Ef, in Sec. 14.11 Finally, Ampere's 
equation, with the current given by Eq. f[2B]) . determines the vector potential in the 
quasiclassical approximation, 

VxVxA{R,t) = / ^Tr{vf{p'f)fsg^{p'f,R,e,t)) . (29) 

Equations f[T5l) - f[22l) and f[28l) - f[29l) constitute a complete set of equations for 
calculating the electromagnetic response of vortices in the quasiclassical limit. For 
high-K superconductors we can simplify the self-consistency calculations to some degree. 
Since quasiparticles couple to the vector potential via ^Vf ■ A, Eq. (!29|) shows that this 
quantity is of order Sne'^Nfv'j/c^ = 1/A^, where A is the magnetic penetration depth. 
Thus, for K = A/^o ^ 1; as in the layered cuprates, the feedback effect of the current 
density on the vector potential is small by factor 

3.2. Linear Response 

For sufficiently weak fields we can calculate the electromagnetic response to linear 
order in the external field. The propagator and the self-energies are separated into 
unperturbed equilibrium parts and terms that are first-order in the perturbation, 

g = go + 5g , A^i=AQ + 5A^i , = ao + 5ai , (30) 
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and similarly for the electromagnetic potentials, A = Aq+6A, <P = 6(P. The equilibrium 
propagators obey the matrix transport equation, 

e - 

[{e+ -Vf Ao)f3-Ao- &o,go]+iVfVgo= 0. (31) 

These equations are supplemented by the self-consistency equations for the mean fields, 
Eqs. ( IT8|) - (fT9|) . the impurity self energy, Eqs. (pT]) - (!22|) . the local charge-neutrality 
condition for the scalar potential, Eq. ( |28l) . Ampere's equation for the vector potential, 
Eq. fl29l) . the equilibrium normalization conditions, 

^o = -^'l, (32) 

and the equilibrium relation between the Keldysh function and equilibrium spectral 
density, 

= tanh [^0^ - ^0^] . (33) 

The first-order correction to the matrix propagator obeys the linearized transport 
equation, 

e ~ 

[{e + -Vf ■ AQ)f3-Ao- ao,5g]^ + iVf -Vdg = [5A ^{ + S&i + 5v , goU , (34) 

with source terms on the right-hand side from both the external field {6v) and 
the internal fields (&4mf, Sai). In addition, the first-order propagator satisfies the 
"orthogonality condition" , 

go'^Sg + Sg'S)go = . (35) 

obtained from linearizing the full normalization condition. Note that the convolution 
product between an equilibrium and a nonequilibrium quantity simplifies after Fourier 
transforming t ^ u: 

Ao^B{e,iu) = Ao{e + uj/2)B{e,iu), B{e,uj) (g) Aq = B{e,uj)Ao{e - uj/2). (36) 

The system of linear equations are supplemented by the equilibrium and first-order 
self-consistency conditions for the order parameter. 



A, {pj,R)=Nj —{V{pf,p'f)f^{p'f,R-e)), (37) 
5A^, {Pf,R;t) = Nj ^iV{pf,p'f)6np'f,R;e,t)), (38) 



and the impurity self-energy, 

&o{Pf,R;e) = niio{pf,pj,R;e) , (39) 

io{Pf,p';, R, e) = u{pj,p'^) + Nj{u{pj,p'^)goip'j, R, ^)Wf,P'f, R, e)) , (40) 

Sai{pf,R; e, t) = n,N}{k{pf,p), R; e) ® 6g{p'f, R, e, t) ® to(P/, P/, R,^)) ■ (41) 

The Keldysh matrix components of the last equation read explicitly 

5af '^(p^, R- e, t) = n,Nf{^'^ipf, p'j, R, e) ® 5^^'^(p}, R; e, t) ® £o^'^(p}, p^, R;e)) , (42) 

6a^iPf, R, e, t) = n,Nf{^{pj,p'^, R, e) ® 5^ (p'^ R, e, t) ® i^{p'f,Pf, R, e)) , (43) 
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'e + uj/2 
2T 

2T 




(44) 



(45) 



In general the diagonal mean fields also contribute to the response. However, we do 
not expect Landau interactions to lead to qualitatively new phenomena for the vortex 
dynamics, so we have neglected these interactions in the following analysis and set 



In what follows we work in a gauge in which the induced electric field, E""^{R; t), is 
obtained from 6^{R] t) and the uniform external electric field, E^^^t), is determined by 
the vector potential SA^(t). For A/,^o ^ 1 we can safely neglect corrections to the vector 
potential due to the induced current. Thus, in the Nambu-Keldysh matrix notation the 
electromagnetic coupling to the quasiparticles is given by 



The validity of linear response theory requires the external perturbation Si) be sufficiently 
small and that the induced vortex motion responds to the external field at the frequency 
set by the external field. At very low frequencies frictional damping of the vortex motion, 
arising from the finite mean free path of quasiparticles scattering from impurities, gives 
rise to a nonlinear regime in the dynamical response of a vortex. This regime is discussed 
extensively in the literature ^0\, and is not subject of our study. However, for sufficiently 
small field strengths the vortex motion is non- stationary over any time interval, although 
it may be regarded as quasi-stationary at low enough frequencies. The non-stationary 
motion of the vortex can be described by linear response theory if 6v 1/r for cj^ 1/r, 
and Si) <^ uj for ujZ^/t. Note that the frequency of the perturbation, uj, is not required 
to be small compared to the gap frequency; it is only restricted to be small compared 
to atomic scale frequencies, e.g u ^ Ef/Ti. 

Self-consistent solutions of Eqs. ( j38l) . (I4T]) and ( l46l) for the self-energies and 
scalar potential are fundamental to obtaining a physically sensible solution for the 
electromagnetic response. The dynamical self-energy corrections are equivalent to 
"vertex corrections" in the Kubo formulation of linear response theory. They are 
particularly important in the context of nonequilibrium phenomena in inhomogeneous 
superconductors. Vertex corrections usually vanish in homogeneous superconductors 
because of translational and rotational symmetries. Inhomogeneous states break these 
symmetries and typically generate non-vanishing vertex corrections. In our case these 
corrections are of vital importance; the self-consistency conditions enforce charge 
conservation. In particular, Eqs. (139|) - (HT|) imply charge conservation in scattering 
processes, whereas (!37j) and ( |38l) imply charge conservation in particle-hole conversion 



A{pj,p'jr) = (i.e. i>inf = 0). As a result the local charge neutrality condition for the 
electro-chemical potential becomes. 




(46) 



Sv = --Vf ■ 5A^(t)f3 + e5^{R] t)l . 



(47) 
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processes; any charge which is lost (gained) in a particle-hole conversion process is 
compensated by a corresponding gain (loss) of condensate charge. It is the coupled 
quasiparticle and condensate dynamics which conserves charge in superconductors. 
Neglecting the dynamics of either component, or using a non-conserving approximation 
for the coupling leads to unphysical results. Self-consistent calculations for the 
local excitation spectrum (spectral density) also provides key information for the 
interpretation of the dynamical response. Because of particle-hole coherence the spectral 
density is sensitive to the phase winding and symmetry of the order parameter, as well 
as material properties such as the transport mean-free path and impurity cross-section. 

In the limit u; ^ the equations above have to be modified, in order to eliminate 
the zero modes associated with the stationary motion of the vortex lattice. However, it is 
possible also to obtain solutions for this limit in quasiclassical theory. For completeness, 
we provide the expression for the transport current density that arises in the low- 
frequency limit, in the flux flow regime [121 IH [13 • 



C 



with 



cell 

"'^^ - lim„ 



de 



-Tr 



e 

2i-v 

c 



■ Aofgjio + -Vf X (V X Ao)f3 



(48) 



5g^ — tanh{e/2T){Sg^ — Sg^) , that is obtained from a systematic 



expansion for small ui. The leading terms are proportional to the vortex velocity vi^. 



4. Nonequilibrium Response of chiral vortices 

The dynamics of the electronic excitations of the vortex core plays a key role in the 
dissipative processes in type II superconductors. Except in the dirty limit, £ -C ^o, 
the response of the core states to an electromagnetic field is generally very different 
from that of normal electrons. It is energetically unfavorable to maintain a charge 
density of the order of an elementary charge over a region with diameter of order of the 
coherence length. Instead, an electrochemical potential is induced which ensures that 
almost no net charge accumulates in the core region. On the other hand, a dipolar- 
like charge distribution develops which generates an internal electric field in the core. 
The internal field varies on the scale of the coherence length. This leads to a nonlocal 
response of the quasiparticles to the total electric field, even when the applied field varies 
on a much longer length scale and can be considered homogeneous. The dynamical 
response of the vortex core includes the collective mode of the inhomogeneous order 
parameter. This mode couples to the electro-chemical potential, 6^, in the vortex 
core region. This potential is generated by the charge dynamics of vortex core states 
and gives rise to internal electric fields which in turn drive the current density and 
the order parameter near the vortex core region. The induced electric fields in the 
core are the same order of magnitude as the external field. The dynamics of the core 
states are strongly coupled to the charge current and collective mode of the order 
parameter. Thus, the determination of the induced order parameter, as well as the 
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spectrum and distribution function for the core states and nonequilibrium impurity 
scattering processes requires dynamical self consistency. Numerous calculations of the 
a.c. response neglect the self-consistent coupling of the collective mode and the spectral 
dynamics, or concentrate on the u ^ limit [TH [291 ED [H!- Presently, quasiclassical 
theory is the only formulation of the theory of nonequilibrium superconductivity capable 
of describing the nonlocal response of the order parameter and quasiparticle dynamics 
in the presence of mesoscopic inhomogeneities and disorder. The numerical solution to 
the self-consistency problem was presented for unpinned s-wave and d-wave vortices in 
Ref. [lOl SH |l2] and for pinned s-wave vortices in Ref. [l3j. In the following we present 
our results for a vortex in a chiral, spin-triplet superconductor. As discussed in our 
accompanying paper [67j, the order parameter we consider is of the form 



Aip,R) = d 



■ J, (Px + ^Pv) ■ J, (Px - "^Pl 

\A+{R)\e''"'^ ^ ^ + \A^{R)\e'P'^^ ^ 



• (49) 



V2 V2 

We have shown, that stable vortex structures can occur for m = —l,p = 1 (singly 
quantized vortex), and for m = —2,p = (doubly quantized vortex with nearly 
homogeneous core). We will in the following discuss for these two vortex structures 
the dynamical charge response, the induced current density, and the magnetic field 
response. Our calculations were performed with impurity scattering included in Born 
scattering limit for a mean free path of i = lO^o, where = hvf/2'KkBTc is the coherence 
length. The Fermi surface parameters are assumed to be isotropic, and the temperature 
was chosen to be T = 0.2 T^. For simplicity we also assumed the high-K limit, where 
the penetration depth is large compared to the coherence length. Our calculations 
of the vortex structure are appropriate to the low-field limit near Hd where vortices 
are well separated from each other. The order parameter, impurity self energy, and 
the equilibrium spectra and current densities were obtained self consistently using the 
Riccati formulation of the quasiclassical transport theory with impurity and pairing 
self-energies [75] . 



Dynamical charge response 

The charge dynamics of layered superconductors has two distinct origins. The c-axis 
dynamics is determined by the Josephson coupling between the conducting planes. Here 
we are concerned with the in-plane electrodynamics associated with the response of the 
order parameter and quasiparticle states bound to the vortex core. We assume strong 
Josephson coupling between different layers, and neglect variations of the response 
between different layers. This requires that the polarization of the electric field be 
in-plane, so that there is no coupling of the in-plane dynamics to the Josephson plasma 
modes. The external electromagnetic field is assumed to be long wavelength compared 
to the size of the vortex core, Aem ^ ^o- In this limit we can assume the a.c. electric 
field to be uniform and described by a vector potential, E^^it) = —^dtA^. We can also 
neglect the response to the a.c. magnetic field in the limit A ^ ^o- In this case the 
spatial variation of the induced electric field occurs mainly within each conducting layer 
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on the scale of the coherence length, ^o- Poisson's equation implies that induced charge 
densities are of order 5$/^q, where 6^ is the induced electrochemical potential in the 
core. This leads to a dynamical charge of order e (A/Ef) in the vortex core. Once the 
electrochemical potential is calculated from Eq. fH6l) we can calculate the charge density 
fluctuations of order {A/EfY from Poisson's equation, 

p^'\R;t) = -^V'6^{R;t). (50) 

Fig. [1] shows the charge distribution for u = 0.2 A which oscillates out of phase 
with the external field. We consider two vortex structures. On the left hand side we 
show results for a singly quantized vortex with quantum numbers m = —l,p = 1, 
and on the right hand side the corresponding results for the doubly quantized vortex 
with quantum numbers m = —2,p = 0. For the singly quantized vortex we observe a 
charge dipole that oscillates with the external frequency. The dipole vector is parallel to 
the external electric field and the charge accumulation concentrated to the vortex core 
area. For the doubly quantized vortex, shown on the right hand side in Fig. [1], a very 
different picture emerges. A dipolar charge distribution accumulates at the interface 
between the dominant and the subdominant order parameter components, oscillating at 
the frequency of the external field. The induced charge which accumulates is of order 
of eA/Ef within a region of order in each conducting layer. As discussed in the 
introduction, this charge is a factor of {6v^Ef) / AP' larger than the static charge of a 
vortex that arises from particle-hole asymmetry [m HSj HH] . 




Figure 1. Out of phase charge response for the case ofanm = — l,p=l vortex 
(left) and for the case of an m = — 2,p — vortex (right) to an a.c. electric field with 
uj — 0.2Z\ polarized in x-direction. Red denotes negative values for the response, blue 
positive. 

Induced current density 

In the static vortex structure a circulating supercurrent is present due to the screening 
of the quantized magnetic field penetrating the vortex. This equilibrium current has 
been calculated in Ref. [67]. For a doubly quantized vortex the circulating current of 
the vortex has a non-trivial structure, with currents in the vortex core that flow counter 
to the circulating supercurrent at large distances from the vortex center. Here, we 
discuss the additional, dynamical part of the current that is induced in the presence of 
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the external a.c. electromagnetic field, and concentrate on the absorptive component 
of the current response, which is in phase with the external field. Results for the a.c. 
component of the current density are shown in Fig. [2] for u = 0.2A. The in-phase 
current response for a singly quantized vortex, shown on the left hand side, shows a 
dipolar pattern, indicating a region of strong absorption {6j\\6E''''^) in the vortex core, 
and emission {6j ■ dE"""^ < 0) in the region roughly perpendicular to the direction of the 
applied field several coherence lengths away from the vortex center. The picture here is 
similar to that for an s-wave vortex in an external a.c. field, discussed in our previous 
work [ini HH HJl |l3j . Energy absorbed in the core is transported away from the vortex 
center by the vortex core excitations in directions predominantly perpendicular to the 
applied field. The net absorption is ultimately determined by inelastic scattering and 
requires integrating the local absorption and emission rate over the vortex array. Note 
that the long-range dipolar component does not contribute to the total dissipation. Far 
from the vortex core the current response is out of phase with the electric field and 
predominantly a non-dissipative supercurrent. 




Figure 2. Absorptive current response for tiie case ofanTO = — l,p=l vortex 
(left) and for the case of an to = — 2,p = vortex (right) to an a.c. electric field with 
Lu ~ 0.2A polarized in x-direction. 

For the doubly quantized vortex, the response of which is shown on the right hand 
side in Fig. [2], we find that the main absorption results from the regions where the 
two time reversed order parameter phases are overlapping. The current response here 
is strongly linked to the order parameter response. Absorptive currents are flowing 
predominantly in the region of the 'domain wall' between the two components. Note 
that the leading contribution in a multipole expansion with respect to the vortex center 
is again the dipolar term, although the pattern here has strong quadrupolar and higher 
order contributions. 

4-3. Magnetic field response 

In the high-K limit we can calculate the a.c. corrections to the local magnetic field in 
the vortex core directly from Eq. ( !29l) . using 

6B{R,t) = V X 6A{R,t). (51) 
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We compare again the singly quantized and the doubly quantized vortex structures 
in Fig. [3l The time- dependent oscillating magnetic field adds to the static, time- 
independent magnetic field that penetrates the vortex in a quantized manner. We find 
that the dynamical magnetic field response is predominantly out of phase, and resembles 
the patterns expected for an oscillation of the magnetic field lines perpendicular to the 
polarization of the external electric field. For the doubly quantized vortex, shown on 
the right hand side in Fig. [3], the magnetic response is concentrated to the region of 
overlap between the two time-reversed order parameter components. 




Figure 3. Out of phase magnetic field response for the case ofanm = —l,p— 1 
vortex (left) and for the case of an m = — 2,p = vortex (right) to an a.c. electric field 
with to = 0.2A polarized in x-direction in the limit of high k. Red denotes negative 
values for the response, blue positive. 

5. Conclusions 

We have discussed the electromagnetic response of a chiral, spin-triplet superconducting 
vortex for two vortex structures, a singly quantized vortex and a doubly quantized 
vortex. The vortex response is nonlocal and largely determined by the response of the 
vortex-core states. We have performed dynamically self consistent calculations, taking 
into account that transitions involving the vortex-core states, and their coupling to the 
collective motion of the condensate are closely related to the dynamical response of 
the order parameter, self energies, induced fields, excitation spectra and distribution 
functions. 

The vortex response roughly consists of an oscillation of the magnetic field lines 
perpendicular to the electric field polarization, and an oscillation of a charge dipole 
parallel to the electric field polarization. Both the charge response and the magnetic 
response are linked to a current fiow pattern that show absorptive response patterns 
in the vortex core. Particularly interesting is the response for the doubly quantized 
vortex. Here we find that the important region is the region where the time-reversed 
order parameter phases overlap, which happens away from the vortex center in a 'domain 
wall' region. 

Based on our results, we see interesting future extensions of our work in the following 
directions: (i) more detailed study of the collective order parameter response for the 
doubly quantized vortex; here in particular the 'optical' mode (oscillation of the two 
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vortex components against each other) is of interest; (ii) to include interactions between 
the vortices by performing calculations on a vortex lattice; and (iii) to study the interplay 
between the charge response and the magnetic response in a vortex lattice. 
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